Vortex states in hole-doped iron-pnictide superconductors 
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Based on a phenomenological model with competing spin-density-wave (SDW) and extended 
s— wave superconductivity, the vortex states in Bai_ a; K a ;Fe2As2 are investigated by solving 
Bogoliubov-de Gennes equations. Our result for the optimally doped compound without induced 
SDW is in qualitative agreement with recent scanning tunneling microscopy experiment. We also 
propose that the main effect of the SDW on the vortex states is to reduce the intensity of the in-gap 
peak in the local density of states and transfer the spectral weight to form additional peaks outside 
the gap. 
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The recent discovery of the iron-pnictide supercon- 
ductors [H, whose parent compounds exhibit long-range 
antifcrromagnetic (AF) or spin-density-wave (SDW) or- 
der similar to the cuprates [2j, provides another promis- 
ing group of materials for studying the interplay be- 
tween magnetism and superconductivity (SC). Espe- 
cially, the hole-doped pnictide superconductors, like 
Bai_ a: K. r Fe2As2 0, have emerged as one of the most im- 
portant systems due to the availability of large homoge- 
neous single crystals. The phase diagram Q for these ma- 
terials indicates that theparent compound, upon cooling 
through Tjv ~ 140K develops a static SDW order. 

By increasing the doping of potassium, the SDW order is 
suppressed and the SC order emerges as the temperature 
(T) falls below T c . The SDW and SC orders coexist in 
the underdoped samples By further increasing the 
potassium concentration to the optimally doped regime, 
the SDW order disappears. These experimental results 
provide compelling evidence for strong competition be- 
tween the SDW and SC orders. 

Another key issue here is the superconducting pairing 
symmetry. Experimental results on the pairing symme- 
try remain highly controversial, leaving the perspectives 
ranging from nodclcss [6, 7] to nodal gap structure 
Although evidence for a nodal gap has been accumulated 
in LaFePO [H and Ba(FeAsi- a; P 2 ;)2 [H3| systems, in the 
K- and Co-doped 122-family of iron pnictides, the ex- 
perimental data points to the existence of isotropic gaps, 
especially in the optimally doped samples Q- Theoret- 
ically it was suggested that the pairing may be estab- 
lished via inter-pocket scattering of electrons between the 
hole pockets (around the T point) and electron pockets 
(around the M point), leading to the so-called extended 
s— wave (ES) pairing symmetry (Ak ~ cosk x +cosk y ) 

In this regard, investigating the vortex states in the 
iron-pnictides can provide useful information on the in- 
terplay between the SDW and SC, as well as the pairing 
symmetry. Recent scanning tunneling microscopy (STM) 



measurement on Bai_ :1 -K.j;Fe2As2 [12| has revealed, for 
the first time in the iron-pnictides, the existence of the 
Andreev bound states inside the vortex core with a sys- 
tematic evolution: a single conductance peak appears at 
a negatively-biased voltage at the vortex center, which 
gradually evolves into two sub-peaks when moving away 
from the center, with a dominant spectral weight at neg- 
ative bias. This negatively-biased conductance peak has 
not been observed in electron-doped Ba(Fei_ a; Co a; )2As2 
[l3| and is beyond current theoretical predictions where a 
peak appears at positive bias in a two-orbital model [l4| 
and at zero bias in a five-orbital model [Hj]. Therefore, 
it is important to develop a sound theory for the vortex 
states in the iron-pnictide superconductors. 

In this work, we adopt a phenomenological model with 
competing SDW and extended s— wave superconductiv- 
ity (ESSC) to study the vortex states in Bai_ 2; K 2 ;Fe2As2 
from the local density of states (LDOS). We show that 
the evolution of the resonance peak in the calculated 
LDOS is in qualitative agreement with STM experiment. 
Moreover, the effect of SDW on the vortex states is also 
discussed, which we predict to be measurable by future 
experiments. 

We begin with a phenomenologically effective two- 
orbital model on a two-dimensional (2D) lattice which 
considers the asymmetry of the As atoms above and be- 
low the Fc layer (lBj . where the onsitc interactions are 
solely responsible for the SDW while the next-nearest- 
neighbor (NNN) intraorbital attraction causes the ESSC. 
The effective mean-field Hamiltonian can be written as 
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(1) 



Here i,j are the site indices, a, j3 = 1,2 are the orbital 
indices, a represents the spin, \x is the chemical poten- 
tial, and rijpa = (Cjp a Cjf3 a ) is the electron density. U 
and Jh are the onsite intraorbital Hubbard repulsion and 
Hund coupling, respectively. Here we have the intcror- 
bital Coulomb interaction U' — U — 2Jh according to 
symmetry [13 • Aij,a0 = S a p^({cj^Cia^) - {cj^c iai .}) 
is the intraorbital spin-singlet ES bond order parame- 
ter, where V is the NNN intraorbital attraction. The 
reason we adopt this model is its ability fiol ] to qualita- 
tively account for the doping evolution of the Fermi sur- 
face and the asymmetry in the SC coherent peaks as ob- 
served by the angle resolved photo-emission spectroscopy 
[HI] and STM [H] experiments on Ba(Fci_ :E Co :r )2As2. 
In the presence of a magnetic field B perpendicular 
to the plane, the hopping integral can be expressed as 

4>/3 = %"/3 ex P[*I^ Ij A ( r ) - dr \i where $ = hc/2e is 
the superconducting flux quantum, and A = (— By, 0,0) 
is the vector potential in the Landau gauge. Following 
Ref. [lal, we have 



t 



h 

2 2 ' 2 t3 

l + (-l) J 4 | l-(-l)^ 
o r 3 H o l 1 

H 




a = P,i = j ±x(y), 
ot = P,i = j ±(x + y), 
a = fi,i =j ± (x - y), 
otj^ P,i= j± (x±y), 
otherwise. 

(2) 

The mean- field Hamiltonian (1) can be diagonalizcd 
by solving self-consistcntly the Bogoliubov-dc Gcnncs 
(BdG) equations: 

H = C^MC, 



(3) 



subject to the self-consistency conditions for the elec- 
tron density and the ES bond order parameter: rijpf = 

J2k=l \Qm-lk\ 2 f(Ek), njpi = 1 — Ylk=l \Qmk\ 2 f{Ek) 
and Ay ij9j8 = -j Y,k=l(Q*mkQnk + Q* n+ x k Qm-\k)f{Ek)- 

Here L = AN x N y , with N x /N y being the number of lat- 
tice sites along x/y direction of the 2D lattice, m = 
Mjy + N y j x ) + 2/3, n = 4(i y + N y i x ) + 2a - 1 and Q 
is a unitary matrix that satisfies (Q^ MQ)k p = SkpEk- 
Here we used i = (i x ,i y ) and j = {j Xl jy), with i x ,j x = 



0,1, 



1 and iy,j y = 0,1,. 



,N,. 



chemical potential \i is determined by the doping con 
centration x through „ X N Y^, 



j0a n iP° 



= 2 - f . The 



ES order parameter at site j is A ■„ = (A 
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A i3 -,^exp[i|- & +3) ' 2 A(r) ■ dr}. The LDOS is given by 
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)/4 where A 



PM = ELlEa [\Qnk?5{uJ-Ek) + \Qn + lk\ 2 5(u + Ek)l 

the supcrcell technique is used to calculate the LDOS. 

In our calculation, the magnitudes of the parame- 
ters are chosen as t-i-i = 1,0.4,-2,0.04 0, U = 3.7, 
V = — 2 and T = 10~ 4 . Magnetic unit cells are in- 
troduced where each unit cell accommodates two su- 
perconducting flux quantum and the linear dimension is 



48 x 24, which is larg er than the coherence 



N x x Ny 

length £ of the iron-pnictides [l3|. Throughout the pa- 
per, the length and energy are measured in units of the 
Fe-Fe distance a and t\, respectively. In the following, 
we focus on two doping concentrations x = 0.4 and 0.3, 
corresponding to the optimally doped and underdoped 
compounds, respectively. 

At x = 0.4, first let us choose Jh = 0.2{7 such that, at 
B = 0, SDW is completely suppressed and the ES order 
parameter A^ is homogeneous in real space. Figures 
0Ja) and ffl^b) show the spatial variations of the reduced 
ES order parameter = A-p{B = 0)| and the 

electron density nj — n jP<? plotted on a 24 x 24 
lattice. The vortex center is located at site (11,12) and no 
SDW is induced. The reduced ES order parameter Aj^ 
vanishes at the vortex center and starts to increase at the 
scale of the coherence length £ to its bulk value, but the 
increase is slower along the 7r/4 and 37r/4 directions with 
respect to the underlying lattice. On the other hand, the 
electron density rij is strongly enhanced at the vortex 
center which is compensated by a depletion of electrons 
around two lattice spacings away from the center, after 
which rij decays also at the scale of £ to its bulk value, 
with no obviously slow variations along the 7r/4 and 37r/4 
directions. The zero-encrgy(ZE) LDOS plotted in Fig. 
[ljc) also peaks at site (11,12) and has the same fourfold 
rotational symmetry (RS) as Aji. In order to reveal the 
spatial variation of LDOS modulated by the vortex, in 
Fig. QJd) we plot the LDOS at four typical positions 
along the black cut in Fig. QJc). As we can see, at the 
vortex center, there is a remarkable negative-energy (NE) 
in-gap peak located at — 0.125A, which is precisely the 
same as observed in Ref. [l2j ]. When moving away from 
the center, the peak will split into two in-gap peaks with a 
dominant spectral weight at negative energy. Finally, the 
LDOS evolves continuously into its bulk feature. The in- 
gap peak and evolution of the LDOS clearly indicate the 
existence of the Andreev bound states inside the vortex 
core, consistent with Ref. fl2l ]. 

In order to study the effect of induced SDW on the vor- 
tex states, we perform the calculation for Jh = 0.23J7. 
Like the Jh = 0-2U case, at B = 0, SDW is completely 



1. The suppressed and the ES order parameter A - g is homo 
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geneous in real space, but the vortex states are funda- 
mentally different from those for Jh = 0.2U and they 
are presented in Fig. [2j As shown in Fig. EJa), the 
vortex center is still located at site (11,12) where the 
reduced ES order parameter Aj^ vanishes, but the size 
of the vortex core is slightly enlarged compared to the 
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FIG. 1: (Color online) Spatial variations of (a) the reduced 
ES order parameter Aja, (b) electron density rij, and (c) 
ZE LDOS map plotted on a 24 x 24 lattice, (d) The LDOS 
at four typical positions along the black cut in (c): at the 
vortex center (11,12); within the vortex core while away from 
the center (10,12); around the edge of a vortex (5,12); and 
far outside a vortex (0,12). The curves in (d) are displaced 
vertically for clarity and the gray dashed line indicates the 
position of zero energy. 



Jh = 0.2C7 case. The induced SDW order parameter de- 
fined as M? = (-iyy(ritf — riji) displayed in Fig. [^b) 
reaches its maximum strength at the vortex center and 
decays at the scale of £ to zero into the superconduct- 
ing region. More interestingly, the SDW order param- 
eter has opposite polarity around two nearest-neighbor 
vortices along the x direction, thus doubling the period 
of the translational symmetry (TS) of the vortex lat- 
tice along this direction. Furthermore, with the induced 
SDW order, the electron density rij is only moderately 
enhanced near the vortex center with a depletion of elec- 
trons around the edge of the vortex core [see Fig. |U[c)]. 
The ZE LDOS in Fig. &ld) also shows a slightly enlarged 
vortex core and the RS inside the core is reduced from 
fourfold to twofold due to the induced SDW order. Fig- 
ures [2je) and[2jf) are the spatial variations of the LDOS 
along the black and white cuts in Fig. d{d), respectively. 
The spectra have only minor differences along the two 
cuts, mainly inside the vortex core and close to lu = 
(indicated by the black arrow in Fig. [H[f)). At the vor- 
tex center, besides a NE in-gap peak at — 0.375A, there 
is a small positive-energy in-gap peak at 0.125 A which 
does not exist for Jh = 0-2U. In addition, the intensity 
of the NE peak is strongly reduced as compared to the 
Jh = 0.2U case and the spectral weight is transferred to 
form additional peaks outside the gap as indicated by the 
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FIG. 2: (Color online) Spatial variations of (a) the reduced 
ES order parameter Afp, (b) staggered magnetization MJ, 
and (c) electron density rij plotted on a 48 x 24 lattice, (d) 
The ZE LDOS map plotted on a 24 x 24 lattice, (e) The 
LDOS at four typical positions along the black cut in (d). (f) 
is similar to (e), but is plotted along the white cut. 

red arrows in Fig. HtL). When moving away from the cen- 
ter, the intensities of all these peaks decrease and finally 
the LDOS evolves into its bulk feature. By comparing 
with the Jh = Q-2U case, we can identify that those two 
in-gap peaks are due to the Andreev bound states while 
the others are due to the induced SDW order inside the 
core. 

At x = 0.3, we choose Jh = 0.32U so that, at B = 0, 
the (7T, 0) SDW coexists with the ESSC. The vortex states 
are plotted in Fig. [31 The reduced ES order parameter 
Ajp [see Fig. |3ja)] shows a y— axis oriented stripe-like 
feature with a modulation period of 8a. The size of the 
vortex core is further enlarged and elongated along the 
y direction. Moreover, the two vortex cores are dragged 
towards each other along the x direction with the vortex 
centers located at sites (15,12) and (32,12), thus also 
doubling the period of the TS of the vortex lattice along 
this direction. The SDW order parameter M| shown in 
Fig. EJb) behaves like nearly uniform stripes oscillating 
with a wavelength of 16a. The vortex core is pinned at 
one of the ridges of SDW stripes where the SDW order 
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FIG. 3: (Color online) (a), (b), (c) and (d) are similar to Figs. 
H», [2tb),[2i;c) and Eld), respectively, (e) The LDOS along 
the black cut in (d). (f) is similar to (e), but is plotted along 
the white cut. 



is stronger than those at other sites. The spatial varia- 
tion of the electron density rij also exhibits a quasi-one- 
dimensional charge stripe behavior with a wavelength 8a, 
exactly half that of the SDW along the x direction [see 
Fig. EIc)]. The one-dimensional stripe-like modulations 
in Aj^. M? and rij already exist at B = 0, which are 
quite similar to the cuprates except for a doubling of the 
period from 4a (for Aj^ and rij) and 8a (for M?) in the 
cuprates [20[ to 8a and 16a in the iron-pnictides. The 
origin of such stripes could be understood in terms of 
the existence of a nesting wave vector oa ~ 0.1257r/a 
connecting the left (right) pieces of the inner and outer 
hole pockets around the T point along the k x direction. 
For proper values of U, Jh and doping, this wave vector 
would modulate AfJ with SDW stripes along the x di- 
rection with period 2n/qA = 16a. The ZE LDOS in Fig. 
EJd) also shows an enlarged vortex core, the doubling of 
the period of the TS of the vortex lattice along the x 
direction and the reduced RS from fourfold to twofold. 
Interestingly, although the two vortex cores are dragged 
towards each other, the ZE LDOS still peaks at sites 
(11, 12) and (36, 12), suggesting that even in the region 
where Aj^ ^ 0, there are ZE states contributing to the 



LDOS. The spatial variations of the LDOS plotted in 
Figs. [He) and Eff) show that at the ZE LDOS peak 
position, there is a NE in-gap peak at — 0.375A, whose 
intensity is further reduced compared to that in Fig. WLX) 
and the intensity decreases when moving away from the 
peak position, indicating that it is due to the Andreev 
bound states. There are also additional peaks outside 
the gap whose positions are similar to those marked by 
the red arrows in Fig. GIL). Their intensities vary dras- 
tically along the black cut in Fig. [3jd) while they barely 
change along the white cut, again suggesting that these 
peaks are due to the SDW order. 

In summary, we have systematically investigated the 
vortex states in Bai_ x K x Fe2As2 with the consideration 
of the interplay between the SDW and ESSC. In the op- 
timally doped compound without induced SDW, there is 
a NE in-gap peak in the LDOS at the vortex center due 
to the Adreev bound states, which splits into two asym- 
metric in-gap peaks when moving away from the center. 
The effect of the induced SDW is mainly to reduce the 
intensity of the NE in-gap peak and transfer the spec- 
tral weight to form additional peaks outside the gap. For 
the underdoped sample where the SDW coexists with the 
ESSC, the vortex cores are dragged towards each other 
along the x direction and the intensity of the NE in-gap 
peak is further reduced. The obtained result at x = 0.4 
without induced SDW is in qualitative agreement with 
experiment and we propose future experiments on the 
near optimally doped and underdoped samples to verify 
the effect of the SDW on the vortex states. On the other 
hand, the disappearance of the Adreev bound states in 
electron-doped Ba(Fei_ 3; Co a: )2As2 may be due to the in- 
duction of strong SDW order in the vortex states, which 
also needs to be verified by future experiments. 
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